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Abstract . Electrostatic charging level of a conducting 
surface in response to injections of electron beams into space 
plasma is investigated by means of one-dimensional Vlasov code. 
Injections of Maxwellian beams into a vacuum shows that the 
surface can charge up to an electric potential <f> s > W b , where W b 
is the average electron beam energy. Since Maxwellian beams have 
extended tails with electrons having energies > W b , it is 
difficult to quantify the charging level in terms of the energies 
of the injected electrons. In order to quantitatively understand 
the charging in excess of W b> simulations were carried out for 
water-bag types of beam with velocity distribution functions 
described by f(V) = A for V min < V < V m?lX and f(V) =0 otherwise, 
where A is a constant making the normalized beam density unity. 

It is found that V max does not directly determine the charging 
level. The pressure distribution in the electron sheath 
determines the electric field distribution near the surface. The 
electric field in turn determines the electrostatic potential of 
the vehicle. The pressure distribution is determined by the beam 
parameters such as the average beam velocity and the velocity 
spread of the beam. 


Introduction 

Electron beam injections from spacecrafts now constitute a 
major activity in space research. Already there are several 
experiments involving rockets ranging in altitude from about 100 
kilometers to about 1500 km (e.g., see Review by Winckler (1980). 
During the STS-3 and Spacelab-2 missions of the shuttle, electron 
beam injections were carrieed out. These space experiments have 
revealed that in response to the injection a host of plasma 
processes are driven (Sasaki et al., 1986; Inan et al., 1984; and 
Shawhan et al., 1984). Low-altitude rocket experiments have 
shown that normally the rocket potential <f> s in response to the 
injection is considerably smaller than the injected electron beam 
energy W b (Winckler, 1980) . This is attributed to the effective 
neutralization of the charges on the vehicle by the return 
current from the ambient ionospheric plasma. At a low 
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ionospheric altitude, where neutral densities are significantly 
high (-10® cm -3 ), beam-plasma discharge provides an additional 
means for charge neutralization. On the other hand, high- 
altitude experiments have demonstrated that the vehicle can 
charge to potentials <f> 9 > W b /e, where e is the electronic charge 
(Managdze, 1983) . The electron beam injections from the shuttle 
showed that when the current-collecting part of the shuttle was 
in the wake and an electron beam was injected, the vehicle 
charged to potentials i> s > W b /e (Sasaki et al., 1986). 

Motivated by the experiments, there are now Several numerical 
simulations on the electron beam injection. Parks et al. (1975) 
used a hydrodynamic approach to study the reflection of 
monoenergetic electron beam injected from a planar body into 
vacuum and showed that the reflection time t r - 2ur b , where 
is the electron-plasma frequency associated with the beam density 
n b . Recently Pritchett and Winglee (1987), Okuda et al., (1987) 
and Okuda and Kan (1987) have used particle-in-cell code to 
investigate the dynamics of the injected beams. Winglee and 
Pritchett (1987) carried out simulations using particle codes 
emphasizing the temporal features of the injected electrons into 
an ambient plasma with density n a « n b . Singh and Hwang (1988) 
carried out simulations using Vlasov codes and dealt with the 
questions of the charging level of the vehicle when the ratio 
n a /n b is varied. They showed that when n a /n b » 1, the return 
current from the background plasma neutralizes the charge on the 
vehicle and the time-average vehicle potential and 0 sa - kT e /e 
where k is the Bolzmann constant and T e is electron temperature. 
On the other hand, when n b > n a , the plasma is not able to 
neutralize the charge and £ sa can appreciably exceed the average 
beam energy. In this case, the electric potential distribution 
near the vehicle is like a thin sheath. In the intermediate 
case, when the vehicle potential lies in the range kT e /e < £ sa < 
W b /e, the beam penetrates into the plasma but much slower than 
the beam velocity. The propagation speed depends on the velocity 
of a triple-charge-layer structure which forms near the 
propagating beam head. 

Despite several simulations, it is not clearly understood why 
vehicles charge to potentials > W b /e (Managdze, 1983; and 
Machlem, 1988) when the ambient plasma is not able to effectively 
neutralize the positive charge on the vehicle. We have 
investigated this issue and we find that the electric field in 
the electron sheath near the vehicle is determined by the 
pressure balance. Since the pressure distribution does not only 
depend on the average beam velocity, but also on the beam 
temperature, density and self-consistent evolution of the plasma 
distribution in the electron sheath, it is difficult to predict 
analytically the dependence of the vehicle potential on the 
injected-beam parameters. Simulations of the injection of 
electron beams with water-bag types of velocity distribution 
functions with a sharp cut-off into a vacuum show that the 
charging level exceeds (m/2e)V„ ax , where V max is the maximum 
velocity above which there are no electrons. This is in contrast 


321 


to the suggestion of Grard and Tunaley (1971) , who assumed that 
for water-bag distributions of photoelectrons emitted from the 
surface, the surface potential z 1/2 (m/e)V m4X . 

The rest of the paper is organized as follows. The numerical 
technique is described in section 2. In section 3, we have 
discussed the injection of Maxwellain and water-bag types of 
electron beams into both vacuum and ambient plasma. The paper is 
concluded in section 4 . 


Numerical Model 

We model the electron beam injection by one-dimensional Vlasov 
simulations, in which the electron beam from a conducting surface 
a t x = 0 is injected into a plasma of extent 0 < X < L as shown 
in Figure 1. The dynamics of the plasma particles and the self- 
consistent electric fields are determined by solving the coupled 
Vlasov and the Poisson equations. The positive charge at X = 0, 
resulting from the injection of the electron beam from this 
surface, is included in determining the electric fields. Any 
charge particles striking the surface is assumed to be lost, but 
their charges are added to the surface charges on the body . The 
surface electric field E x (x = 0) is determined by the net surface 
charge density. The plasma particles which exit the boundary at 
X = L are reflected back into the system, simulating a uniform 
plasma. When the beam electrons begin to reach the surface at X 
= L, the simulation is stopped. At X = L, we use the Dirichlet 
boundary condition *(X = L) =0, which is found to be good as 
long as the perturbations created by the electron beams do not 
reach this boundary. 

In the simualations described here we use L = 10 \ d , where \ d 
is the Debye length with a reference plasma density n 0 and 
temperature T e . We have used the electron to ion mass ratio for 
H + . The numerical grids in X and V x space are as follows: aX = 
X., aV x for electrons is 0.25 V. and for ions 0.05V t , where V^ is 
the electron thermal velocity with the temperature T e . The time 
step to advance the solutions is At = 0.1 up 0 , where u po is the 
electron— plasma frequency with the density n^ . Further details 
about the simulation technique can be found in the work of Singh 
and Schunk (1984). , a 

We hav£ used the foll 9 wing normalizations: distance X - X/\ d , 

velocity V = V/V t , time t = tw p? , current density J = J/J t , J* = 
en 0 V t , density n = n/n 0 , potential $ - e^/T e , electric field E - 
E/E , E 0 = T e /e\ d . In the case of electron beam injection into 
a vacuum, n 0 and T e refer to the injected beam. When the beam 
injection occurs into an ambient plasma n 0 and T e refer to 
the ambient plasma. 


Numerical Results 

We begin this section with the discussion on the injection of 
electron beams into a vacuum. Such an exercise throws a great 
deal of light on the causes for the vehicle charging considerably 
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in excess of W b /e, where W b is the average electron beam energy, 
W b = l/2m # V b . The results on the injection into an ambient 
plasma is described in the following section. 


Injection Into Vacuum 

iMaxwellian Beams . The temporal evolution of the spacecraft 
potential for several velocities of Maxwellian beams are shown in 
Figure 2, in which V b is the beam velocity normalized to its 
thermal velocity V t . During the early stage (t < 3u pb ) the 
potential increases at a fast rate and then it settles down at a 
quasi-steady value depending on the average beam velocity. It is 
seen that the quasi-steady value of 0 8 slowly increases. This is 
attributed to the fast electrons in the tail of the Maxwellian 
beam which continually escapes. The analytical calculation of 
Grard and Tunaley (1971) show an infinitely large potential when 
the ambient plasma is completely absent. The vehicle potential 
saturates at a time approximately given by 


This time is somewhat longer than the beam reversal time t r : 

2 wl b calculated by Parks et al (1975), using a hydrodynamic 
treatment for the beam propagation. However, it should be noted 
that the t r is the time when the beam velocity V b -» 0 in the 
retarding potential distribution which t 8 is the time when the 
reflected beam electrons reach the surface at x = 0 and 
effectively neutralize the further increase in the positive 
charge on it. 

Figure 2 shows that the surface vehicle potential increases 
with the average beam velocity V b . The dependence of $ 8 on V bA at 
t = 3iip b is plotted in Figure 3, which also shows the plot of W b 
versus A V b . Comparing the two curves in this figure we find that 
jL » w b . Intuitively it can be argued that in a Maxwellian A beam 
yith a finite temperature there are electrons at velocities V > 

V b and therefore attains a value for which such fast electrons 
are also confined by the developing electric fields. But in a 
Maxwellian distribution there is no unique maximum velocity which 
can uniquely determine the maximum possible value of the surface 
potential £ 8 . 

Water-Baa Beams . In order to quantitatively understand the 
dependence of $ t on the energy of the injected electrons, we 
carried out simulations with water-bag type of distribution 
functions defined by 
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where A is chosen so that the bean density is unity, namely, 
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The water-bag distribution functions have the attractive property 
that they have a sharp cut-off at V = V Bax with no electrons at V 
> v max . Thus, simulations with such distributions can possibly 
show* the dependence of ^ s on the maximum electron energy in the 
beam, which is not a well-defined quantity for a Maxwellian beam. 

The solid curves in Figure 4 show the temporal evolution of 
the surface potential f 9 for two water-bag beams with the same 
average beam velocity V b = 6 A but V Bax and A V B j n for the two beams 
are different. Beam f 1 has V max = 12 and V Bin = 0, while for 
beam #2 V Bax * 7 and V B l n =■ 5. ^Beam #1 is warm while beam #2 is 
relatively cold. For beam #1, $ 8 = 170 while for beam #2 p 9 = 

75. These values far exceed the charging level determined by the 
maximum kinetic energies of the electrons in the two beams. 

These energies are 


* „I 

- V,. x = 72 , 
2 

beam #1 

(4) 

«..x - 2-1-5, 

beam #2 

(5) 


Grard and Tunaley (1971) have carried out analytical 
calculations on the charging of a plane surface in response to 
photoelectron emissions. They have considered a water-bag 
distribution for the emitted electrons. They have suggested that 
the charging level is determined by the maximum electron energy 
in the velocity distribution function. Our simulations show that 
this is not true as the surface charges to a potential greatly in 
excess of W Bax (Figure 4) . This authors used continuity and 
energy conservation equations 

nV = n 0 V 0 (6) 

1 1 

— mV 2 - e.<p = — mV 2 - e£ 0 (7) 

2 ’ 2 0 

and the Poisson equation to show that the surface electric field 
is given by 
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where f 8 (V) is the distribution function on the surface and N s is 
the electron density there. In equations (6) and (7) n and V are 
the density and velocity where electric potential is 0 and 
similarly n„ and V„ are the density and velocity where ^ = 0 O . 

We now compare these anlaytic results with our results 
obtained from the numerical simulations. In our normalized units 
equation (8) can be written as 
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2N 8 



f 8 (V)dV 


( 9 ) 


where f 8 (V) is one-sided distribution function and it is 
normalized to take into account both outgoing and incoming 
particles. We note that the integral in (9) is the effective 
electron temperature near the surface, if the average drift 
velocity is zero. 

Figures 5a and 5b show the electron velocity distribution 
functions near the surface for beam #1 and #2 respectively. The 
effective electron temperature associated with the distribution 
functions is given by 
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where u is the average drift velocity associated with the 
distributions. We find that near the surface for both the A beams 
T eff is about 25 (see Figure 7), and the electron density N s = 2 
(see Figure 6). Substituting these values in (9), we obtain E ? = 
10, which is remarkably close to the electric fields obtained m 
the numerical simulations as shown in Figure 4 by the dotted 
curves. We find that despite this fair agreement on the surface 
electric fields obtained by theory and simulations, the surface 
potentials (^ 8 ) obtained from simulations greatly exceed those 
predicted by the theory (Grard and Tunaley, 1971) . We now show 
that this difference is caused primarily by the pressure 
distribution in the electron sheath near the surface. We find 
that (7) is not valid throughout the electron sheath because it 
does not include the thermal energy. A A 

The spatial distributions of density A Nj(X) , effective 
temperature T eff (X) and electric field E(X) in the simulations 
for beam #1 and #2 are shown in Figures 6, 7 and 8, respectively. 
Using these distributions we now examine the relative 
contributions of the terms in the momentum balance equation, 
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where u is the average drift velocity and P is the pressure given 
by 

P = N T e (12) 

If u = 0, the electric field distribution is given by 
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n ax 
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In Figure 8, the solid and dashed 9 urves show the numerically 
obtained spatial distributions of E for the two beams. The curve 
with the triangles is the electric field obtained from (13) for 
the beam #1. The pressure in (13) is obtained by the density and 
temperature shown in Figures 6 and 7, respectively. It is seen 
from Figure 8 that the electric field is quite accurately 
deteripingd £y the pressure force. This shows that the kinetic 
term u 3u/3x is negligibly small in the momentum equation. This 
is expected from the nature of the distribution function in 
Figure 5a for beam #1. A A A 

The curve with circles (o-o) gives the plot of N' 1 3P/3X for 
beam #2. It is seen that it slightly underestimates the electric 
field shown by the dashed curve. This is accounted by the fact 
that for the beam #2 the distribution function has distinct peaks 
(Figure 5b) , which is not described well by an effective 
temperature . 

The comparison of the electric fields directly obtained from 
the simulations with those derived from the pressure force 
clearly shows that it is the pressure force in the electron 
sheath which determines the electric field distribution and hence 
the surface potential of the vehicle. Comparing the electric 
field curves for the beam #1 and #2 in Figure 8, we see that the 
beam with the large V max has relatively large electric fields 
extending to much greater distances than the beam with the 
smaller V max . Thus, there is a dependence of charging on V max , 
but it is not directly determined by the maximum electron energy 
1/2 m. V„ ax . The dependence is determined by the pressure 
distribution in the electron sheath. It is difficult to predict 
the pressure distribution because of the highly non-linear nature 
of the problem. 


Effects of Ambient Plasma on Charging 

The ambient plasma reduces the charging level by providing a 
return electron current, which partly neutralizes the positive 
charge on the surface. Figures 9a and 9b show the temporal 
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evolution of the surface potential for a Maxwellian beam with 
beam velocity V b = 6 injected into ambient plasmas with densities 
n ■ 0.1 n b and 0.2 n b . The beam and ambient plasma temperatures 
are assumed to be the same. The effect of the ambient plasma on 
the surface potential is found to be twofold? it causes 
oscillations in the surface potential and the time average value 
of 4> decreases with increasing ambient plasma density. Without 
any ambient plasma, the surface potential attained a quasi-steady 
value of about ? 8 = 70 (see Figure 3). For n a /n b = 0.1 and 0.2, 
its time-average values reduce to 50 and 45, respectively. 

However, it is important to note that for such small ambient 
plasma densities, th$ charging level exceeds considerably the 
average beam energy W b = 18, and the beam does not propagate into 

the plasma. . . . , 

The oscillations seen in Figure 9 in the surface potential are 
found to have time periods determined by the ambient plasma 
frequency. For example, the time periods seen in Figures 9a and 
9b are-r = 20 vZl and 14w nK , respectively, and they are given 
2* ->J(n /n ) ' 

b a pb 

We have found that for sufficiently small beam densities when 
4> < W b , the beam begins to penetrate into the plasma (Singh and 

Hwang, 1988) . However, the propagation speed critically depends 
on the time-average surface potential 0 sa . When - (kT/e) « 

W b for n a » n b , the beam propagates into the plasma with the 
injection velocity Singh and Hwang (1988). As ^ 8 increases with 
decreasing value of the relative ambient plasma density (n a /n b ) , 
the propagation velocity decreases. For surface potentials in 
the range (kT/e) « * s < W b , Singh and Hwang (1988) found that 
the beam head penetrates into the plasma with well-defined 
laminar potential structure near the beam head. The velocity of 
the potential structure is determined by the nonlinear plasma 
processes through which it evolves. 

Conclusion 

We have investigated in this paper the charging level of a 
conducting surface when an electron beam is injected from it. 
Injections into both vacuum and ambient plasma are considered. 
When a Maxwellian beam is injected into vacuum it is found that 
the surface charges to a potential much greater than the average 
beam energy. The dependence of the charging in excess of the 
average beam energy is investigated by considering beams with 
water— bag types of velocity distribution functions with distinct 
maximum velocity V max such that no electron velocities V > V max . 
It is found that the electric field distribution in the electron 
sheath near the surface is determined by the pressure 
distribution. Thus, the surface potential is determined not only 
by the V nax , but by all the beam parameters such as the density 
n b , average beam velocity and the velocity spread of the beam. 

The parameters determine the pressure distribution in the 
electron sheath. Since the effective temperature and the density 
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distributions in the sheath evolve self-consistently through the 
nonlinear dynamics of the electrons, it is difficult to predict 

the charging level. , , . 

In our one-dimensional model, we have not included return 
currents coming from directions other than that for the beam 
injection. For the omnidirectional return current, the surface 
potential will in general tend to be smaller than that for the 
one-dimensional case. However, space experiments have shown that 
vehicles do charge to potentials comparable to, or in excess of, 
the average electron beam energy (Sasaki, et al. 1986, Managdz, 
1983, Machlem, 1988). Therefore, the physical processes seen in 
the one-dimensional simulations are relevant to space 
experiments. 
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right hand boundary the electric potential #(x~d) « o 
while the left hand boundary (x - 0) is Kept floating 
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A ^ 

Figure 3. Surface potential as a function of beam velocity V b 
for the beam§ injected into vacuum. Normalized beam 
energy (1/2 V 2 b ) as a function of V b is also shown. 



"v 


Figure 4. Temporal evolutions of surface potential 0, and 

electric field E 8 are sh 9 wn for water-bag be^ms #1 and 
#2. For both the beams V b = 6; for bfam #1 V Bln = 0, 
V Bax 58 12 and for *> eam #2 , V Bin » 5, V B4X = 7. 
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(b) Beam #2 
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(a) Beam #1 
X = 0+ 



Figure 5. Electron velocity distribution functions near the 

surface X = 0 + at tw pb * 50. (a) beam #1, (b) beam #2. 

It is worth mentioning that the distribution functions 
shown are not exactly at x = 0, but they are at x = 0* 
as they are averaged over 2 grid spaces in front of the 
surface at x = 0. The distributions show that the 
injected water-bag beams have been modified by the 
intense electric fields near the surface. 



Figure 6. Density distribution in the electron sheath for the 
water-bag beams #1 and #2. 








Figure 7. Effective temperature distributions for the beam #1 and 
#2 in the electron sheath. 



Figure 8. Electric field distributions. Beam #1: solid curve 
gives E fyom numerical simulations and the triangles 
give the E field fr 9 m pressure balance. Beam #2: 
brokep curve gives E from simulation while the circles 
give E field from pressure balance. 
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Spacecraft Pot. vs. Time 


n a /rib = 0.1 


n«/nh - 0.2 


Figure 9. Temporal evolution of surface potential 0, when 

Maxwellian beams are injected into ambient plasmas: 

(a) ambient plasma density n a =0.1 n b , (b) n a = 0.2 n 
beam velocity V b = 6. 






